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Resonant inelastic x-ray scattering (RIXS) lias become a powerful tool in the study of the elec- 
tronic structure of condensed matter. Although the linewidths of many RIXS features are narrow, 
the experimental broadening can often hamper the identification of spectral features. Here, we 
show that the Maximum Entropy technique can successfully be applied in the deconvolution of 
RIXS spectra, improving the interpretation of the loss features without a severe increase in the 
noise ratio. 



I. INTRODUCTION 

Over the last decade or so, the technique of resonant 
inelastic soft x-ray scattering (RIXS) has developed into 
a formidable tool in the study of the electronic structure 
of the solid state [1, 2], nanomaterials [3], and even liquids 
and gases [4, 5]. The strengths of RIXS using soft x-rays 
owe to its bulk sensitivity, atomic (and even orbital) se- 
lectivity, ability to sample finite momentum transfer and 
narrow linewidth. Indeed, the recent development of sub- 
100 meV resolution grating spectrometers [6] has opened 
up the new possibility of using soft x-rays to study low- 
energy collective excitations such as magnons or orbitons, 
and even to track their dispersion in momentum [7, 8] . 

In the RIXS process, a core electron is excited into an 
unoccupied valence orbital, creating a core hole in the 
intermediate state which then rapidly (and coherently) 
decays to the final state via x-ray emission (see, for ex- 
ample, Refs. [1, 2]). The incident x-ray is tuned near the 
threshold of a particular core electron excitation, grant- 
ing RIXS its site and orbital selectivity. For example, 
for transition metal L-edges the process is of the form, 
2/ 3d" -^ 2p5 3d"+i -^ 2p6 3d"*, in which the * denotes 
an excited state. The final state may be an excited con- 
figuration (e.g. dd* crystal-field transition) or a collective 
excitation, and the energy difference between the incident 
and emitted x-rays represents the excitation energy. 

Typically however, for soft x-ray RIXS where charge- 
transfer and/or crystal-field excitations are of interest, 
the combined energy resolution of the incident photons 
and RIXS spectrometer amounts to 0.5 - 1.0 eV near the 
O K-edge (^ 520 eV). On the other hand, the typical 
energy of transition metal dd* crystal field excitations is 
of the order of 1 - 4 eV, and their separation can be close 
to the limit of the resolving power of moderate-resolution 
instruments (for example, low-energy dd* transitions in 
VO2 [9, 10]). Moreover, even in high- resolution measure- 
ments, different low-energy collective excitations can lie 
in close proximity relative to the resolution function; for 
example, bi- and single magnon peaks as well as phonon 
contributions within 400 meV in La2Cu04 [11]. There 
is therefore a sensitive trade-off between statistical pre- 



cision and resolution for these kinds of measurements, 
which are perfectly poised to take advantage of the ben- 
efits of a reliable deconvolution procedure. 

The technique of Maximum Entropy (MaxEnt) owes 
its origins to the study of communication theory intro- 
duced by Shannon [12], in which the proposed measure 
of information content, S, had the same form as the ther- 
modynamic entropy. 
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in which pij is the number of counts in a pixel (i, j) 
and k is an arbitrary constant. The basic idea of Max- 
Ent is relatively simple: one maximizes the information 
content, 5", of the processed (deconvoluted) probability 
distribution subject to it being consistent with the mea- 
sured data. This consistency test is achieved through a 
X^ comparison between the processed distribution, con- 
voluted with the (known) resolution function, and the 
measured distribution. The outcome of this process is 
a distribution that is 'most likely' to have been respon- 
sible for the measured data, given the properties of the 
resolution function. 

The algorithm employed here is the Cambridge Algo- 
rithm [13] and has been successfully applied to many 
fields of data analysis, e.g. positron annihilation [14, 15], 
image analysis, astrophysics and extended x-ray absorp- 
tion fine structure data [16]. The solution of the algo- 
rithm is iteratively updated in three 'search directions' 
di.j,m with coefficients a„: 
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The search directions used in this algorithm are VS* (to 
maximize the entropy), Vx^ (to minimize the x^), and a 
third direction involving higher derivatives. The solution, 
p[ ■ , and coefficients a„ are updated at each iteration 
until convergence in the solution is achieved. The par- 
ticular benefit of the MaxEnt procedure over more tra- 
ditional Fourier transform-based deconvolution methods 
is the substantially improved signal-to-noise ratio (which 



will be discussed in more detail in Section V), even when 
faced with sparse (or missing) data [16]. 

In the following, all measurements were performed at 
the AXIS endstation of beamline 7.0.1 of the Advanced 
Light Source, Lawrence Berkeley National Laboratory, 
employing a Nordgren-type grating spectrometer [17]. 
Spectra were recorded on a 2D multi-channel plate detec- 
tor. Correction for the curvature introduced by the opti- 
cal components of the spectrometer was achieved by fit- 
ting the peak position of several well-characterized sharp 
emission features across the non-dispersive axis of the de- 
tector. This correction was then used for all subsequent 
spectra. The MaxEnt procedure was applied to the raw 
2D spectra before curvature correction. This was found 
to produce more favorable results compared with pro- 
cessing corrected, integrated ID spectra, a finding that is 
expected owing to the greater information content of the 
2D images. The MaxEnt deconvolution process involves 
convolving the processed data with the known resolution 
function (the 'broadening' function), for which there are 
two components in RIXS measurements: i) the energy 
resolution of the incident photons, and ii) the spectrom- 
eter resolution. The second of these, predominantly aris- 
ing from the finite source size, is much more easily dealt 
with: it is approximately Gaussian with respect to the 
wavelength of the emitted photons. An accurate knowl- 
edge of the incident photon resolution is much more chal- 
lenging since its impact on emission spectra is non-trivial 
and depends in part on the specific excitations involved 
in the vicinity of the incident photon energy (i.e. the ab- 
sorption spectrum); no attempt has been made to remove 
this component from the experimental data. 



II. DECONVOLUTION OF REFERENCE 
SPECTRA 

As an initial diagnosis of the performance of the Max- 
Ent procedure, x-ray emission spectra of the L3^2-edge 
of Zn were obtained at various different spectrometer slit 
widths, corresponding to different spectrometer resolu- 
tions. These spectra represent transitions of the form 
3d -^ 2p with the spin-orbit split L3 edge at a lower 
energy than the L2 emission line. Since we are dealing 
with emission features, excited well above the absorption 
threshold, the incident photon energy does not contribute 
to the overall resolution of these measurements, and they 
therefore provide a robust test of the MaxEnt procedure 
in removing the spectrometer part of the resolution func- 
tion. In Fig. 1, Zn spectra recorded (in second order [18]) 
with spectrometer resolutions of 0.28 eV and 0.78 eV at 
FWHM are shown, alongside the results of the MaxEnt 
deconvolution of the 0.78 eV spectrum. A broadening 
function of 0.66 eV (85% of the spectrometer resolution) 
was used in the MaxEnt deconvolution routine. It is clear 
from Fig. 1 that the width of the deconvoluted spectrum 
is much narrower than the raw spectrum, and indeed ap- 
proaches the width of the narrower 0.28 eV spectrum. 



r = 0.78 eV, raw data 
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FIG. 1: (Color online) Comparison between raw (second or- 
der) Zn x-ray emission spectra recorded at 0.78 eV spectrom- 
eter resolution and the results of the MaxEnt deconvolution 
procedure. For comparison, a spectrum recorded at 0.28 eV 
spectrometer resolution is also shown, a) Integrated spectra, 
b) raw 2D spectrum for F = 0.78 eV, c) MaxEnt deconvoluted 
2D spectrum for F = 0.78 eV, and d) raw 2D spectrum for 
F = 0.28 eV. The raw 2D spectra (before correcting for the 
curvature of the image) are shown on a logarithmic false color 
scale. 



indicating that a large portion of the instrument resolu- 
tion has been removed from the spectrum. In fact, this 
narrowing of the line width of the emission lines is already 
visible in the raw 2D spectra shown in Fig. Ib-d. The 
FWHM of the L3 edge (including the natural width of 
this feature) is 1.29 eV for the raw spectrum, compared 
with 1.06 eV after deconvolution. For comparison, the 
FWHM of the same feature in the narrower raw spectrum 
is 1.00 cV. Moreover, there are no additional artefacts in- 
troduced in the deconvolution procedure - the MaxEnt 
spectrum closely resembles the 0.28 eV spectrum. The 



behavior of the deconvolution with varying broadening 
functions was also investigated, and found to be very 
stable for functions of FWHM < 90% of the total res- 
olution function. Above this, some artificial sharpening 
close to the emission peaks was observed. Finally, it is 
worth noting that the signal-to-noise ratio is not signifi- 
cantly decreased in the deconvoluted spectra (a subject 
to which we will return in Section V). 



III. DECONVOLUTION OF CRYSTAL FIELD 
EXCITATIONS 

In order to test its performance in resolving close spec- 
tral features, the MaxEnt deconvolution procedure was 
applied to Co L-edge RIXS data of C03V2O8, which is a 
kagome staircase compound consisting of Co^+Og octa- 
hedra separated by V^''"04 tetrahedra [19]. Spectra were 
recorded in second order with a spectrometer resolution 
of 0.82 eV and an incident photon resolution of 0.4 eV, 
amounting to a total resolution of approximately 0.91 eV. 
For the MaxEnt deconvolution procedure, a broadening 
function of 0.82 eV was used, equivalent to the spectrom- 
eter resolution. The same procedure was applied to Co 
L-edge emission reference spectra in order to check that 
no artefacts were introduced in the deconvolution. 

The raw RIXS data are shown on a loss energy scale 
in Fig. 2a for several different incident photon energies 
between 777 eV (spectrum a) and 782 eV (spectrum g), 
spanning the Co L^-edge absorption feature. The peak at 
eV represents elastically scattered light. There are clear 
loss features present in all spectra below —4 eV, above 
which a weak and broad charge-transfer peak emerges 
centered about —6.7 eV. In order to accurately locate 
these peaks, a linear combination of Voigt functions was 
fitted to each spectrum, and the average center of each 
feature was determined. Altogether, three distinct loss 
features can be seen in the spectra: (i) at —0.83 eV 
clearly visible in spectra (d) and (f), (ii) at —2.08 eV 
in spectra (c-f), and at —3.37 eV in spectrum (d). These 
spectra, and the energies of these features, are similar to 
RIXS L-edge measurements of other Co^"*" compounds, 
for example CoO [20]. 

The MaxEnt deconvoluted spectra are shown in 
Fig. 2b, and show obvious improvement in the linewidth 
of the features. In these spectra, the loss features previ- 
ously identified are now much more clearly visible. For 
example, features (i) and (ii) are directly visible in all 
spectra, and (iii) is visible in all but the lowest energy 
spectrum (a). It is emphasized that no attempt to re- 
determine the energy location of these features with the 
new information provided by the MaxEnt algorithm has 
been made, yet they agree within ~ 0.1 eV in all spec- 
tra. It is noted that the noise level of the MaxEnt spec- 
tra is slightly higher than for the raw data, as might be 
expected in any deconvolution procedure. Nevertheless, 
focussing on features that persist across several incident 
photon energies allows one to be sure that the origin of 
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FIG. 2: Co L-edge RIXS data of Co3V208. (a) Raw spec- 
tra, and (b) MaxEnt deconvoluted spectra, in which the raw 
spectra are reproduced in light grey. Solid lines are a guide 
to the eye and represent the data with a binomial smoothing 
of 0.15 eV. Vertical lines show the elastic line at eV and 
loss features identified in the raw data. The symmetry of the 
excitations are labelled in (a). 



a feature is intrinsic to the system under study, and not 
an artefact of the increased noise fioor. The strength of 
this approach is anticipated to lie when spectral features 
are difficult to identify in raw data, as for example in 
spectrum (c), in which the —0.83 eV feature is hard to 
separate from the elastic peak. 



IV. DECONVOLUTION OF ELASTIC PEAKS 

As an additional test, the MaxEnt procedure has been 
applied to V La-edge RIXS data of NdVOa at room tem- 
perature, a system in which the occupation of the V 3d 
orbitals becomes ordered at low temperature [21]. Spec- 
tra were recorded with a spectrometer and incident pho- 
ton resolution of 0.36 eV, yielding a total energy resolu- 
tion of around 0.51 eV. Seven spectra were obtained ap- 
proximately equally-spaced across the V La-edge absorp- 
tion feature (and are shown in Fig. 3), and each spectrum 
was treated to the MaxEnt procedure with a broadening 
function of 0.36 eV. Again, the same procedure was ap- 
plied to a Zn L-edge emission reference spectrum to en- 
sure no artefacts were introduced in the deconvolution. 
The relative intensity of RIXS spectra at the V L-edge 
is weak, and so the spectra here have been summed to- 
gether to yield 'averaged' loss spectra, shown in Fig. 3d 
for both the raw and MaxEnt deconvoluted spectra. The 
idea behind this procedure is that fluorescent features, 
that are dispersive in loss energy, contribute a broad and 
weak background, whereas loss features will reinforce in 
the summed spectra. Note that there will be a slight ad- 
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and provide a weak background to the summed spec- 
tra presented in the manuscript. In the raw summed 
spectra (Fig. 3d), two principle features are present: the 
elastic peak at O eV and a broader feature at around 
—2 eV. Shoulders either side of this second peak indicate 
the presence of other spectral features. Once the Max- 
Ent procedure is applied, however, the location of these 
weaker features becomes clear, at —2.62 eV and —1.09 eV 
either side of the —1.85 eV peak, and represent crystal 
field dd* excitations. These features are also clearly iden- 
tifiable in the individual deconvoluted spectra shown in 
Fig. 3c, particularly the lowest energy spectrum. Fur- 
thermore, however, a relatively weak shoulder is evident 
close to the elastic peak at —0.42 eV, and may represent 
an orbital excitation in the form of a bi-orbiton, previ- 
ously observed for YVO3 at 0.4 eV [22]. Its presence 
can also be inferred in the raw spectra from the slight 
asymmetry of the elastic peak in this spectrum. Finally, 
an additional feature at —3.74 eV, not visible in the raw 
spectrum, becomes clear after the MaxEnt deconvolu- 
tion. This peak is harder to identify in the individual 
spectra, but is most prominant in the highest energy 
MaxEnt spectrum in Fig. 3c. In order to check its pres- 
ence at more than one excitation energy (and ensure that 
its origin is not the fluorescent part of the spectrum), the 
highest energy MaxEnt spectrum has been removed from 
the summation in the dotted line of Fig. 3d, and indeed, 
this peak still persists. However, a higher energy feature 
(at around —4.5 eV) disappears in this process; although 
this may be a loss feature of the spectrum, further work 
is needed to confirm its origin. 
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FIG. 3: (Color online) V L-edge RIXS data of NdVOs. 
(a) The absorption spectrum and excitation energies used for 
RIXS. Raw RIXS data are shown on (b) an emission energy, 
and (c) a loss energy scale, for which the individual MaxEnt 
spectra are also shown, (d) Summed RIXS spectra before 
and after MaxEnt deconvolution. The symmetry of the ex- 
citations are labelled in (d); the dotted line represents the 
MaxEnt summed data with one spectrum missing (see text). 



ditional broadening of the features due to uncertainty in 
the initial photon energy. 

The individual raw spectra are shown in Fig. 3b. The 
first spectrum is on the onset of the V La-edge, in a 
location most suitable for exploring loss features. At 
higher excitation energies, V 3d fluorescence begins to 
contribute more strongly in the spectra, but is mostly 
concentrated beyond 2.5 eV from the elastic peak. There- 
fore, the subsequent interpretation of features in the 
summed spectra below this energy is less influenced by 
the presence of fluorescence. Moreover, the contribu- 
tions from fluorescence do not reinforce across spectra. 



The strong elastic peak in these data mean that we 
can directly measure the instrument resolution, since the 
elastic peak is a (5-function in the limit of an infinitessi- 
nially small total resolution function. In order to avoid 
complications with low-energy (e.g. phonon) excitations, 
all the following fits to the elastic peak have been con- 
strained predominantly to its high-energy side. For the 
raw spectra, the FWHM of the elastic peak is 0.75 eV, 
slightly broader than expected, presumably in part due 
to the additional processing involved in summing these 
spectra (but possibly also due to some uncertainty in 
the slit width). For the MaxEnt spectra, however, the 
FWHM is 0.50 eV. Although at first glance this is not 
as narrow as one might expect, it should be remembered 
that this represents a combination of the incident pho- 
ton resolution and the spectrometer resolution, and no 
attempt has been made to remove the incident photon 
resolution function. In this respect, the narrowing of the 
elastic peak behaves well, and corresponds to a narrowing 
of the 'effective' spectrometer resolution (by subtracting 
in quadrature the incident photon resolution, 0.36 eV, 
from the measured resolution) from 0.66 eV to 0.35 eV. 



V. NOISE PROPAGATION 

The propagation of noise through the MaxEnt proce- 
dure is a complex problem; indeed, with only a slightly 
different setup of the problem, MaxEnt can be used to 
'de-noise' data or reconstruct missing information from 
sparse data [16, 23]. In order to quantify the propagation 
of noise through the MaxEnt procedure, we have simu- 
lated a series (of M — 100) of noisy RIXS spectra, Pi{E). 
These have then been passed through the MaxEnt decon- 
volution, and the variance in the resulting deconvoluted 
spectra analyzed. 
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where P{E) is the mean of the M simulated spectra. 
The test spectrum, T{E), was chosen to approximate 
spectrum (d) in Fig. 2: a linear combination of Gaus- 
sian functions of 0.4 eV FWHM (similar to the beamline 
resolution used in the measurements) centered at eV, 
—0.83 eV, —2.08 eV and —3.37 eV were added to a Gaus- 
sian function of 4 eV FWHM centered at —6.73 eV to ap- 
proximate the elastic, dd* and CT features respectively. 
This spectrum was then convoluted with a Gaussian func- 
tion of 0.82 eV FWHM, approximating the spectrometer 
resolution, and scaled to contain 1000 counts in the peak 
data channel. The M simulated spectra were then cre- 
ated by adding randomly generated noise following a nor- 
mal distribution with a{E) — y/T{E), and are shown in 
Fig. 4a. These test spectra, Pi{E), were each deconvo- 
luted with the same parameters as used for the C03V2O8 
data shown in Fig. 2. 

The resulting deconvoluted spectra are shown in 
Fig. 4b. Comparison between the deconvoluted spectra 
and the initial test spectrum, T{E) [shown by the dashed 
line in Fig. 4b], are very encouraging, and provide a di- 
rect visualization of the power of the MaxEnt algorithm. 
In order to ensure that the mean and variance are mean- 
ingful quantities, the inset shows the distribution of the 
data about the mean, in units of the standard deviation, 
for data between —9 and 1 eV (to restrict the contribu- 
tion to the finite signal region of the spectra). Although 
there is a weak positive skewness in the distribution of 
the deconvoluted spectra (emphasizing the complexity of 
the MaxEnt noise problem), it is sufficiently close to a 
normal distribution that the mean and variance are still 
useful indicators of the distribution. This is important 
since it allows us to attach meaningful (statistical) errors 
to the deconvoluted spectra. 

It is evident from Figs. 4a, b that, as expected, the noise 
level is slightly higher for the deconvoluted spectra (for 
example, compare the error bars for the peak at eV in 
Fig. 4a and at —2 eV in Fig. 4b, for which the intensity 
is similar). For Poisson statistics, relevant in counting 
problems, the variance of a datum scales with its ex- 
pected value (for x > 10), as shown in Fig. 4c in which the 
variance, a'^, of the input data is plotted against its mean. 
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FIG. 4: (Color online) The results of deconvoluting M = 100 
simulated noisy spectra with the MaxEnt procedure; (a) Input 
spectra, and (b) deconvoluted spectra. The individual spec- 
tra, Pi{E), are shown in light gray, with the mean spectra, 
P{E), in blue. The error bars shown represent the standard 
deviation, a[P{E)], of the M spectra for each energy. The 
inset shows the distribution of each simulated spectrum from 
its mean, normalized to its standard deviation. In (b), the 
initial test spectrum, T{E), is shown by the dashed line, (c) 
Comparison between the mean of the input spectra and their 
variance, a^. (d) Comparison between the mean of the de- 
convoluted spectra and a^'^'''. For (c) and (d) the solid line 
represents P = a^ and P = a^'^^ respectively. 



P . However, this is not the case for the deconvoluted 
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FIG. 5: (Color online) The (a) raw, and (b) MaxEnt decon- 
voluted RIXS spectra shown in Fig. 2, reproduced here with 
error bars. In (a), the statistical error is shown {a = vN), 
whereas in (b) the empirical error a = A'^^'^'^^ is used (see 
text). 



spectra, indicating the propagation of errors through the 
MaxEnt process is non-hnear with respect to the pixel 
intensity. QuaUtatively, the relationship is super linear, 
meaning that pixels of high intensity (in the deconvoluted 
spectra) are more sensitive to statistical noise than those 
of lower intensity. Such behavior is connected with the 
tendency of the deconvolution to move counts from low 
intensity regions of the input spectrum towards higher 
intensity regions. Empirically, we find that the variance 
and mean are connected by P ^ cr^-^^ for the range of P 
investigated here (approximately 50 - 2000), as shown in 
Fig. 4d, which corresponds to an approximate doubling 
of the noise ratio for a typical spectrum. It is empha- 
sized that this analysis only reflects the propagation of 
statistical noise, and does not account for systematic er- 
rors that may be present in the process itself. The result 
of applying these empirical errors to the C03V2O8 RIXS 
data is shown in Fig. 5 for some representative spectra 
before and after the MaxEnt procedure. In each case. 



the visual scatter of the data points is consistent with 
the magnitude of the error bars, and the loss features 
that were previously identified are well above the noise 
level. Moreover, the apparent structure at high energies 
in the deconvoluted spectrum (g) is of the order of the 
noise, and is due to the poorer statistics recorded for this 
spectrum. 



VI. CONCLUSIONS 

In summary, the MaxEnt deconvolution procedure has 
been successfully applied to soft x-ray RIXS spectra. The 
deconvoluted spectra show a marked improvement in the 
resolution of spectral features without introduction of 
artefacts associated with the process or excessive increase 
in the noise ratio, allowing for greater confidence in sep- 
arating and identifying loss features, such as crystal-field 
dd* transitions or low-energy collective excitations. For 
example, a very recent application of MaxEnt has helped 
to clarify the RIXS features of Lai_2;Sr2;Mn03 [24]. De- 
tailed analysis of the propagation of noise through the 
MaxEnt procedure has been presented, and the noise ra- 
tio of deconvoluted spectra has been found to approxi- 
mately double for the features of typical spectra. The 
process is quite general, and is not limited to the soft x- 
ray regime or 2D data, and is expected to perform equally 
well with, for example, x-ray absorption spectra or high- 
resolution RIXS. It is anticipated that the strength of 
the procedure is when spectral features are difficult to 
identify in raw data (but whose presence may already be 
inferred from that data, albeit indirectly). 
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